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A REVIEW OF HIGH-SPEED, CONVECTIVE, HEAT-TRANSFER 
COMPUTATION METHODS 

Michael E. Tauber 
NASA Ames Research Center 
Moffett Field, California, U.S.A. 

Abstract 

The objective of the paper is to provide useful engineering formulations and to instill a 
modest degree of physical understanding of the phenomena governing convective aerodynamic 
heating at high flight speeds. Some physical insight is not only essential to the application of the 
information presented here, but also to the effective use of computer codes which may be available 
to the reader. The paper begins with a discussion of cold- wall, laminar boundary layer heating. A 
brief presentation of the complex boundary layer transition phenomenon follows. Next, cold- wall 
turbulent boundary layer heating is discussed. This topic is followed by a brief coverage of sepa- 
rated flow-region and shock-interaction heating. A review of heat protection methods follows, 
including the influence of mass addition on laminar and turbulent boundary layers. The paper con- 
cludes with a discussion of finite-difference computer codes and a comparison of some results 
from these codes. An extensive list of references is also provided from sources such as the various 
AIAA journals and NASA reports which are available in the open literature. 

Nomenclature 

B mass-addition parameter (Eq. 70) 

C mass concentration 

Ci constant in Eq. 37 

Cf skin-friction coefficient 

Cp specific heat at constant pressure 

D binary diffusion coefficient 

F diffusion function (Eq. 6) 

f similarity stream function (Eq. 13) 

g total enthalpy ratio, H/He 

H total enthalpy, h + (u 2 + w 2 )/2 
H cavity depth (Eq. 65) 

h static enthalpy 

h ’ reference enthalpy (Eq. 53) 

j exponent (j = 0 for yawed, infinite cylinder; j = 1 for axisymmetric body) 

This paper is declared a work of the U.S. Government and is not subject to copyright protection in 
the United States. 
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k coefficient of thermal conductivity 

k roughness element height (Eq. 50) 

L cavity length (Eq. 65) 

Le Lewis number, pCpD/kf 

M Mach number 

m mass addition rate 

in molecular weight of gas 

Nu Nusselt number, q w c p x/k(h aw - h w ) 

Pr Prandtl number, (iCp/k 

p pressure 

Pp peak pressure (Eq. 68) 

p s stagnation pressure 

q w heat transfer rate to the body surface 

r radius of body of revolution 

r n nose radius of hemisphere 

r enthalpy recovery factor 

Re Reynolds number, pux/p 

St Stanton number, cjw/peUeOw - h w ) 

t static enthalpy normalized by total enthalpy at boundary layer edge, h/H e 

u,v,w velocity components in the x,y,z directions, respectively 
Voo ffee-stream velocity 

x,y,z chordwise, normal, and spanwise coordinates 
a angle of attack 

P pressure gradient parameter (Eq. 21) 

T heat-transfer function ratio (Eq. 43) 

Y specific heat ratio 
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8 cone half angle, wedge angle, or flat-plate angle with respect to free stream 

e inverse density ratio across shock, poo/p e 

0 boundary layer momentum thickness 

A sweep, or yaw, angle of wing leading edge or cylinder 

|i coefficient of viscosity 

M-0 reference value of coefficient of viscosity 

^,T) similarity variables (Eq. (9-10)) 

p density of gas 

(p density and viscosity product ratio (Eq. 17) 

\\f stream function (Eq. 11) 

\\r ratio of heat transfer at surface with mass addition to value without mass addition (cold- 
wall value) (Eq. 69) 

C spanwise velocity ratio, w/w e 

Subscripts 

A,M atoms and molecules in binary mixture 
aw adiabatic wall (no heat transfer at wall) 
bt beginning of boundary layer transition 

c cone 

cyl cylinder 

e value at boundary layer edge 

f chemically frozen value 

FP flat plate 

s stagnation point value 

SL stagnation line 

w wall, or local body surface, value 
°o free-stream conditions 
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Introduction 

The original impetus for developing means of predicting the aerodynamic heating during 
high-speed atmospheric flight was the advent of the long-range missile. Only a modest extension 
of technology was required from the missile atmospheric entry velocity of about 6 km/sec to the 
7.5-km/sec entry speed of vehicles returning from satellite orbit. It is the fortuitous behavior of air 
molecules at high temperatures that made the early formulations of aerodynamic heating applicable 
to a speed of about 9 km/sec, thereby covering both the missile and satellite entry flight regime. 

The pioneering papers of Lees 1 and Fay and Riddell, 2 published in 1956 and 1958, 
respectively, provided the basic physical and mathematical modeling for hypervelocity convective 
aerodynamic heating. However, the early work 1 * 2 was hampered by a lack of information of the 
high-temperature thermodynamic and transport properties of air. Succeeding computations by 
Beckwith and Cohen^ and Cohen, 4 published in 1961, used the high-temperature thermodynamic 
properties of Moeckel and Weston 5 and the transport properties of Hansen, 6 both published in 
1958. The state of the art in convective heat transfer at hypersonic speeds in the early 1960s was 
well presented by Dorrance. 7 Possibly the best current source of information for the working 
engineer is the chapter on forced convection in external flows by Rubesin et al. 8 

In a high-velocity flow, the physical mechanism of heat transfer differs from the molecular 
heat-conduction process occurring at lower speeds. At sufficiently high speeds, the air that has 
passed through the strong shock wave surrounding the front of the body becomes so hot that the 
molecules dissociate into atoms. When the atoms diffuse through the boundary layer into the cooler 
region near the wall, the heat of recombination that is released as molecules form can contribute 
significantly to the heating of the wall. Therefore, at high speeds convective heating consists of 
both conduction and diffusion (due to atomic recombination) of energy through the boundary 
layer. The ratio of energy transported by diffusion to that by conduction is known as the Lewis 
number and is an important parameter in high-speed, laminar, boundary layer calculations. The 
chemical reactions in the laminar boundary layer are frequently frozen at altitudes in excess of 
about 45 km according to Fay. 9 Therefore, the heating can be potentially reduced 2 if the wall is 
made of a noncatalytic, glassy, or ceramic material which inhibits atom recombination. For exam- 
ple, Rakich et al. 1 ® have shown that the Shuttle orbiter tiles are noncatalytic. Metals and many 
ablative heat-shield materials that are widely used on missiles are catalytic; therefore, finite-rate 
catalytic effects will not be considered in the following discussions. Wall catalysis effects are 
included in the finite-difference boundary layer and surface material interactions computer codes 
which will be mentioned later. 

The objective of the present review is to provide useful engineering formulations and to 
instill a modest degree of physical understanding of the phenomena governing convective aerody- 
namic heating. Some physical insight is not only essential to the application of the information pre- 
sented here, but also to the effective use of computer codes which may be available to the reader. 
(In many codes, the physics are obscured by the mathematical procedure required to solve the dif- 
ferential equations.) An extensive list of references is also provided from sources such as the vari- 
ous AIAA journals and NASA reports which are available in the open literature. 

The paper begins with a discussion of cold- wall, laminar boundary layer heating. (The term 
"cold-wall" means that no mass addition occurs at the body surface and that the influence of the 
wall affects the heating only through the temperature of the wall.) A brief presentation of the com- 
plex boundary layer transition phenomenon follows. Next, cold-wall turbulent boundary layer 
heating is discussed. This topic is followed by a brief coverage of separated flow-region heating. A 
review of heat protection methods follows, including the influence of mass addition on laminar and 
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turbulent boundary layers. The paper concludes with a discussion of comparisons and some results 
from finite-difference computer codes. 


Within the text of the paper, the MKS metric system of units is used. However, most of the 
figures that are shown originated in other papers and are not replotted if English units were 
originally used. Therefore, conversion factors are presented in the Appendix. 

Laminar Heat Transfer — Cold Wall 

In this section, analytic solutions are derived for the laminar boundary layer heating of sur- 
faces with and without pressure gradients. Air will be treated as a real gas in equilibrium and con- 
tinuum flow is assumed. However, the changes in heating that occur in low-density flows will be 
discussed. 

Boundary Laver Equations 


The boundary layer equations are written in the usual form, which assumes that the thick- 
ness of the layer is small compared to the body's radius of curvature and that centrifugal forces can 
be neglected. Following the historic development of Refs. 1 and 2, the equations are further sim- 
plified by assuming a binary mixture of "air molecules" (instead of nitrogen and oxygen) and "air 
atoms in thermodynamic equilibrium. (The binary gas mixture assumption greatly simplified the 
computation and was initially justified because the final results agreed well with experiments. Sub- 
sequently, it was theoretically validated by Moss. 11 ) The expressions for the conservation of mass, 
momentum in the x and y directions, and energy are 4 



The coordinate systems are shown in figure 1. For a three-dimensional body of revolution, 
j - 1 and the velocity component w = 0. For a two-dimensional, infinite, yawed cylinder' 
j — 0, w = w e = constant and 9w/9z = 0 in the boundary layer. 

For a binary gas mixture in thermodynamic equilibrium, the continuity equation does not 
have to be written for the individual species because the species concentrations are uniquely deter- 
mined by the pressure and enthalpy, or any other two state properties. In the energy equation the 
diffusion function, F, is defined as 


F(p,h) = (Le-l)(h A -h M )| 


9h 


4>=const. 


( 6 ) 
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The boundary conditions for the conservation relations, Eq. (1-5), are: 
at y = 0 


u = v = w = 0 
H = H w (x) 


( 7 ) 


and at the boundary layer edge 


u = u e (x) and w = w e = const. 
H = He = const. 


( 8 ) 


The momentum and energy equations, Eq. (3) and Eq. (5), respectively, are now trans- 
formed to modified Lees similarity coordinates 1 ’ 2 by letting 


£00 = “T J PwlJ-wUe^-j^ dx 


and 


ri(x,y) = 


^0 0 

Ue(r/LV 


mV 2(5 - 5) ' 0 




p dy 


(9) 


( 10 ) 


where po is a reference value of the coefficient of viscosity, and ^ is a function of x, as yet 
undetermined. To satisfy continuity, Eq. (1), the stream function is introduced 



Combining Eq. (10) and Eq. (11a) gives 



Now a new dependent variable can be defined as 


f(x,h) 


¥ 


no V 2(5 - 5) 


where the velocity profile is given by 


Df u_ 
3r| “ Ue 


(11a) 

(Hb) 


( 12 ) 


(13) 


( 14 ) 


Writing the nondimensional spanwise velocity and enthalpy ratio, respectively, as 
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g(^ ,T| ) = 


(15) 

(16) 


and the ratio of the products of density and viscosity as 


<P 


= _£±L 

Pe4e 


(17) 


permits transforming Eq. (2, 3, and 5) into the similarity coordinate systems. Now the momen- 
tum and energy equations become, respectively, 




a 2 f 

arj2 




sr( ,p $■} 


+ 2(E-E)fiLJ5L ns , 

K M an d<]di ~ si a n 2 1 <18) 




d_ 

9ti 


.ftf 0 + F) l^ 


f _\ 

dj[ 


1 - 




J 




£ I& 1 - Prf + F) b - (i J + (1 - « sr« 2 >]} 




where the pressure gradient parameter, P, is defined as 


(19) 


( 20 ) 


p =2i§ 

Ue te d£ 


(21) 


and the ratio of static to total enthalpy is 


t_ He~ g_(1_ts) ^ 2 " (ts ~ te) ^J 


( 22 ) 


For the body of revolution at 
constant and given by 


a - 0, t s — 1 and £ — 0. On the infinite yawed cylinder, t s is 


w 2 V 2 sin 2 A 

ts=1 -2He= 1 - ! k- 


( 23 ) 
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which is the static enthalpy ratio at the boundary layer edge along the stagnation line. The boundary 
conditions in the similarity coordinate system for Eq. (18-20) are at q = 0 (wall) 


f(o,o = |^-(o,^) = C(o^) = o 

(24a) 

g(0»4) = gw(^) = tiv(^) 

(24b) 

as h — > °° (boundary layer edge) 


|^(~4) = C(“4) = g<~4) = i 

(25) 


Similar Solutions 


Partial differential equations were traditionally difficult to solve numerically. However, 
when written in the local similarity coordinates, the equations can be readily solved numerically. 
Exact solutions can also be found for a (limited) number of conditions where the flow external to 
the boundary layer and the wall temperature became independent of the streamline coordinate 
^(x). Practical examples of such exact similar solutions are the flow at a three-dimensional stagna- 
tion point, at the stagnation line on a two-dimensional infinite cylinder, and on an inclined sharp 
flat plate or a sharp cone having a constant wall temperature. For the preceding conditions, the 
partial differential equations become ordinary differential equations, which are much easier to 
solve. 


Most practical problems of interest, however, do not have exact solutions. A widely used 
approximation is to assume that such flows can be treated as if they were locally similar. The pro- 
cedure consists of assuming that the derivatives with respect to ^(x) vary relatively slowly and can 
be neglected, while the local values are used for the terms which depend on £(x) through the 
external flow or the wall conditions. Again, a set of ordinary differential equations results which 
can be solved for the required boundary layer profiles. The solutions thus found are then applied at 
each value of £(x) in the nonsimilar flow by using the approximate values of _t s and the local val- 
ues of te, g w » and (3. Under the local similarity conditions, the value of d^/dE, becomes arbitrary 
and may be set equal to zero. 3 The momentum and energy equations, Eq. (18-20), can now be 
written as ordinary differential equations, where primes denote differentiation with respect to q, to 
yield 


(<pf ')' + fT = p 


JL 


^(l + F)g- 


.< r > j - {^y - 't - h\ 

(q>C’y + fC* = 0 

+ fg* = fe(l - Prf + F) [(ts - te)(f)’ 2 + (1 - t s )(C 2 )’] 


(26) 

(27) 

(28) 


The boundary conditions are 


f(0) = f ’(0) = C(0) = 0 
g(0) = g w 


(29a) 

(29b) 
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while the parameters (5, t s , te, and g w are evaluated locally along the body. The similar boundary 
layer equations can be numerically integrated in a straightforward manner for high-speed flight 
conditions if the high-temperature thermodynamic and transport properties are known. 

Thermodynamic and Transport Properties 


In 1958, the first set of high-temperature thermodynamic^ and transport^ properties for air 
became widely available and were extensively used for the following two decades. Although more 
precisely calculated values for the transport properties were published in the early 1960s by Peng 
and Pindroh 12 and Yos, 12 the reports were not widely circulated. In Ref. 14-16, the transport 
properties from Ref. 6, 12, and 13 were compared and it was shown that Hansen's coefficients 
of viscosity and thermal conductivity were too low at temperatures above 2000 K. While such 
inaccuracies should not be ignored, it was shown by Howe and Sheaffer 17 that varying the coef- 
ficient of thermal conductivity by a factor of 10 resulted in a 40% increase in the stagnation point 
heating rate at a speed of about 18 km/sec. 


For use with computer codes, the analytic approximations for high-temperature air proper- 
ties developed by Worbs and Bolster 18 are useful. When fast computers are available, the program 
of Gordon and McBride 19 can be used as a subroutine and gives thermodynamic properties and 
species concentrations that compare closely with those of Ref. 5. 

Heat Tra nsfer Rate at a Stagnation Point or Line 

The general expression for the heat transfer rate to the wall can be written 4 

qw = ^ (1+F » ) (|r), (30a) 

or, in the similar coordinate system 


q »=i^< I + F w)g'w(e^ d £j‘ 


(30b) 


The dimensionless heat transfer parameter consisting of Nusselt number and Reynolds number 
evaluated using static wall conditions, can be written as J 


where 


Nu w g' w 

(Re w )0- 5 gaw — gw 


(1+F) 


L x dUe^P-5 
|3te Ue dx J 


(31) 


gaw - te + r(l — te) 

and r is the recovery factor. Two exact solutions to Eq. (26-28) consisting of the axisymmetric 
stagnation point and the stagnation line on a yawed, infinite cylinder will be presented next The 
equations were numerically integrated using Ref. 5 and 6. The calculations were performed to 
temperatures up to the onset of ionization, which corresponds to a flight speed of almost 9 km/sec. 

A unit Lewis number was used so that F = 0; this assumption can introduce errors of up to 8% in 
the heating rate. 4 K 

, For *' “symmetric stagnation point j = 1, p = 0.5, t s = t e = 1 and due/dx = Ue /x, 
so that Eq. (31) gives & 
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Nu w _ V2 g w 
(Re w ) 0 - 5 ” 1 - gw 

Correlation of the numerical calculations yields 4 


(1+F) 


^ Uw - 0 767Pr^‘ 4 f P e ^ e 

(Re w )0-5 " °* 767Pr w f [p w ^ w j $ 


0.43 


which can be compared to the expression from Ref. 2 of 


NU y 

(Rev 


^^=0.768Pr°’ 4 f-£^ 

; w ) 0 - 5 w f (.PwPw J s 

Using Eq. (32) yields the heat transfer rate 


0.40 


qw s = (H e - hw)s(pePe)°- 43 (PwPw) 0 - 07 


( 32 ) 


(33) 


(34) 


For the yawed, infinite cylinder j = 0, P = 1, t s = t e < 1, and due/dx = Ue/x so that 

(35) 


Nu w g'w /1 , p \ 
(Re w )°- 5 gaw — gw w 


For a recovery factor, r, of 0.85 and 


PePe 

PwPw 


< 1 


Ref. 4 gives the expression 


Nu w 

(Re w ) 6 - 5 


= 0.57Pr 


0-4 

-w f 


- \0.45 

PcPe ^ 

PwPwJ, 


(36a) 


and for 


■ P#e } 
PwPw J 


> 1 


Nu w 

(Re w )°- 5 



(36b) 


where Eq. (36b) is applicable for large yaw angles and highly cooled walls. Note that for 

(_£etPL| = i 

lPwPwJ s 
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the ratio of Eq. (32) to Eq. (36a) disagrees by only 5% from the theoretical value of 2' 0 - 5 , which 
relates the stagnation heating rates on an axisymmetric body and a cylinder. 


The usefulness of the correlation equations are further improved by writing the relations 
directly in terms of flight condition. Examples of such correlations are shown in figure 2 for the 
ammar heat transfer parameter for the stagnation point on a hemisphere, an unswept circular 
cylinder, sharp cones, and a sharp flat plate. (A wall temperature of 300 K was used for the cal- 
culations illustrated in figure 2 to facilitate comparison with shock tube data.) Note that the maxi- 
mum difference between the four hemispherical stagnation point calculations 2 - 2 *}' 22 j s a 5 out \$% 
and can partly be attributed to using different high-temperature transport properties. The sharp 
cone values 2 - 3 - 24 differ by 8% at low speeds and much less at high speeds. 

The calculations shown as straight lines in figure 2 can be used to derive simple analytic 
approximations for the heat transfer rate. Noting that 



(Re w )0-5 ^ 


(37) 


where the exponent "a" varies from 0.3 to 0.4, "b" varies from 0.10 to 0.20, and Ci is a constant 
and substituting the expressions for Nusselt and Reynolds number, yields the axisymmetric stag- 
nation point heat transfer rate 6 



- h 


w 


V -3 p 0-5 



(38) 


where p s is the stagnation pressure. The velocity gradient at a hemispherical stagnation point can 
be evaluated using Newtonian theory, 2 which gives 


where 



r n 


2(Ps - Poo)' 


0.5 


Ps 


(39a) 


Ps - PooV^I - + p 0 

and e = p«./p s . At high speeds e « 1, and 



(39b) 


An empirical expression for e 
10 3 m/sec to 10 4 m/sec, is 


which is valid within about ±10% in the speed range from 


e = 11.4V 


-0.57 


OO 


(39c) 
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For wall temperatures of less than about 4000 K, the Prandtl number variation is small and the 
coefficient of viscosity can be approximated by 


Hw~ 


t 0.71 


If the exponent, b, in Eq. (37) is set equal to 0.145, then the bracket in Eq. (38) becomes inde- 
pendent of wall temperature. Using a value for the exponent, a, of 0.3575, and if the flight speed 
is high enough so that hoo « V^/2, one can now write Eq. (38) in the simple and useful form 


q Ws » 1.83(1(H) 



(40) 


where the nose radius, r n , is in m, the free-stream density is in kg/m 2 , and the flight velocity is in 
m/sec. The correlation equation for air given by Marvin and Deiwert 21 is nearly the same as 
Eq. (40). (Stagnation point heating correlations in gases other than air are given in Refs. 21 
and 25.) The analogous expression for the swept, infinite cylinder is 

q Wcyl -1.29(l(H)^J' 5 (l-0.18sin2A)V^l - A ^ (41) 


where 


h aw - hoc + 0.5V^(1 - 0.18 sin 2 A) 

and A is the sweepback angle. In writing Eq. (41) it has been implicitly assumed that the stagna- 
tion point velocity gradients are identical for the sphere and cylinder. According to Ref. 8, the 
Newtonian value, Eq. (39), is accurate for the circular cylinder, but 7.7% too low for the sphere. 
If the velocity gradient value of Ref. 8 were used in Eq. (40), the constant would increase from 
1.83 to 1.90. 

A comparison with experimental data is shown in figure 3 for the four hemispherical stag- 
nation point correlations of figure 2. Up to a speed of 14 km/sec, all four correlations, from 
Ref. 2 and 20-22, lie within the data spread, which is about ±25%. In addition to the shock tube 
data, five measurements based on observations of the onset of melting of projectiles fired in a 
ballistic range by Compton 26 are also shown. The ballistic range data confirm the shock tube 
measurements and theories of Refs. 20 and 21 in the low ionization velocity regime. At speeds 
above 14 km/sec where ionization effects are very strong, Hoshizaki's formulation 20 appears to 
give the best agreement with the data of Rose and Stankevics 27 The strong ionization regime is 
also characterized by intense radiative heating; fortunately, the speeds are much higher than those 
of concern here for missiles. 

Newtonian theory, Eq. (39), can be used to predict the velocity gradients reasonably well 
on circular cylinders and hemispheres. However, for blunter configurations more precise numeri- 
cal methods, or experimentally derived values, must be used. In figure 4, a comparison between 
experimentally determined velocity gradients and values calculated using Newtonian theory shows 
the limitations of the Newtonian approximation for determining stagnation point velocity gradients 
on blunt bodies. 



Blunt-Bodv Heating-Rate Distributions 
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The distribution of the heat-transfer rate around a blunt-nosed axisymmetric body, or a 
yawed cylinder, can be determined using the local boundary layer similarity approximation. Fol- 
lowing Ref. 4, the ratio of the local heat transfer to the stagnation value is written as 


Qw _ 


qw s 


0.5 


^PwM-w du e ^ / ( PwM-w du e ^ 

t Pte dxj / [ pte dxjj 


(42) 


where the velocity gradient parameter, p, is 


P = 


2t s (due/dx) 


PwPwUgteCr/L) 2 )^ - Q 


. _ I PwPwUe [' 


PwPwUe I r 2 dx 


(43) 


and 


r = 


(\ + F 


w 


Pr 


w 


>,/ 


1 + F w 
Pr 


W 


g w s 


The gradient of the enthalpy ratio at the wall, g' w , is determined by solving the locally similar 
boundary layer Eqs. (26-28). The full solution of Eqs. (42) and (43) requires iterating, because 
F is inside the integral in Eq. (43). However, it is consistent with the local similarity assumption 
to set T = 1 in Eq. (43), thus eliminating the iteration. The result is essentially the same as that 
of Kemp et al 28 and is shown in figure 5 for a hemisphere and figure 6 for a circular cylinder. 
Also illustrated in figures 5 and 6 is the solution of Ref. 1, which results when T = 1 in both 
Eq. (42) and (43). For the axisymmetric body, numerical solutions were correlated to yield 


Qw 

qwc 


2pwPw | ■JJ”’ 


-|0.5 


lK)L(2£)0-5 1 + 0.096(Btp)°- 5 

PwPwUer 1.068 


(44) 


for a Lewis number of one and g w « 1. 

Also shown in figures 5 and 6 are heating solutions following a cosine variation of the 
body angle measured from the stagnation point. Although it is approximate, the cosine variation of 
heating has a theoretical basis (Ref. 1) at least to the body sonic point (near 45°) and appears to 
give reasonable results to about 70° beyond the stagnation point. 

The local similarity solution procedure has also been applied to a blunt-faced axisymmetric 
body with a small comer radius. The results of using the methods of Refs. 1 and 28 to calculate 
the heating distribution are illustrated in figure 7. Again, the method of Ref. 28 is seen to give 
somewhat better results than that of Ref. 1 , although there is much scatter in the data. ~ 
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Cone and Rat Plate Heating Rates 


The laminar boundary layer heating of bodies without pressure gradients 
((3 = due/dx = 0) such as sharp cones and flat plates at angle of attack will now be considered. 
The expression for the heating rate on a sharp cone in air is, from Ref. 23, after noting that the 
total enthalpy term should be replaced by the recovery enthalpy, 


<lwc = 


0.018 





(45) 


where p e is in atm, u e is in m/sec, x is in m, and h is in m 2 /sec 2 . At high speeds, where 
2 

hoo « V J2, and using the Newtonian flow approximation 



and the transport properties of Ref. 15, Eq. (45) becomes 



4.03(10-5) 


^Qoo cos 8 c ^ q - 5 v 3.2 

x J °° 


sin 


8r 1 


KA 

haw^ 


W 

m 2 


(46) 


where 


h aw - hoc + 0.40V 2 

“ oo 


and 8 C is the cone half angle. In practice, all cones have some tip bluntness. Usually, the nose 
blunting will reduce the heating on the cone surface aft of the nose. In an analysis made for a per- 
fect gas of y = 1.4, Vallerani 2 * showed that the local heating is reduced at a flight Mach number 
of 5, or greater, for 5° < 8 C < 25°, but increased for larger half-cone angles. 

The sharp wedge, or inclined flat plate, heating rates can be found from Eq. (45) by using 
the Mangier transformation which is 


q WFP - Qw c /(3)® - 5 (47) 

However, the Newtonian approximation for p e u e is not as accurate for a wedge as for a cone. 
After adjusting the Newtonian values, an approximation for the inclined flat plate which is equiva- 
lent to Eq. (46) is 


4wfp“ 2.42(10-5)^ 


COS 


SfpY 


.5 


r 3.2 


W 


v"- z sin Spp 1 - iF- (48) 


m* 


where 8pp is the wedge angle, or angle of attack of the plate. 

Wing Leading-Edge Heating 

The heating along a cylindrical leading edge of a finite-length wing can be approximately 
calculated using the expression of Rubesin 12 
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<3wle 


Q + Q 

M W C yl ^WFP 


sin 2 A 


N0.5 


cos a 


(49) 


where qw cy j is given by Eq. (41) and qpp by Eq. (48). For highly swept leading edges, 
Eq. (49) is limited to small angles of attack, since the stagnation line moves from the cylindrical 
leading edge onto the lower surface of the wing as the angle of attack increases. 


Low-Densitv Flow Heating 


The peak hypersonic aerodynamic heating experienced by missiles or manned vehicles 
always occurs in the continuum flow regime. However, sometimes vehicles have small nose or 
wing leading-edge radii, and these regions of the vehicle can experience low-density flow 
phenomena which alter the heating. For example, if a vehicle flying at high altitude has a small 
nose radius, the Reynolds number on the nose ceases to be large. In that case, the basic assump- 
tions of the Prandtl boundary layer theory are violated as the boundary layer becomes thick and slip 
can occur at the body surface. In addition, shock waves can become thick, or diffuse, and occupy 
- a significant fraction of the shock layer. The low-density flow field and heating are discussed by 
Shorenstein and Probstein 3 ^ and Shorenstein 32 for flat plates and cones. Axisymmetric stagnation 
point heating measurements are presented and compared with theory by Boylan. 33 One parameter 
that characterizes low-density flows is the Knudsen number, which is the ratio of molecular mean 
free path to body size. The increase in stagnation point heat transfer coefficient with Knudsen 
number and altitude was calculated by Moss et al. 34 and is shown in figure 8. For a nose radius 
of 2.54 cm, the heat-transfer coefficient, predicted by the direct simulation Monte Carlo (DSMC) 
method, increases tenfold in the altitude range of 60-100 km. Of course, the heat- transfer rate 
decreases by two orders of magnitude because of the atmospheric density decrease with altitude. 
From figure 8, it is also evident that neglecting slip at the wall can lead to a 130% overprediction 
of the heating at a Knudsen number of 0.1. (The results from DSMC methods have been verified 
by comparison with experimental data by Nomura. 35 ) 

This concludes the discussion of cold-wall laminar heat transfer. Next, boundary layer 
transition will be briefly covered. 


Boundary Laver Transition 

The transition of the boundary layer from laminar to turbulent flow remains the most com- 
plex problem in fluid mechanics. At hypersonic speeds, turbulent boundary layer heating can be 
several times greater than laminar heating. Therefore, it is essential that some form of reliable 
means of predicting transition be available to avoid the penalties that result from overly conserva- 
tive design. Because transition is influenced by many factors, the engineer must rely on empirical 
relations derived from test data. However, none of the ground test facilities can simulate most of 
the parameters of interest; in fact, the operating characteristics of many test facilities have been 
found to strongly influence the data. In the following sections, some examples of test data will be 
discussed and a few correlation charts and formulas from the open literature will be presented. 

It has long been known that among the important parameters influencing transition are the 
boundary layer edge Reynolds number and Mach number. The results of plotting measurements of 
the Reynolds number calculated at the beginning of transition against Mach number are illustrated 
in figure 9 for flow on cones. At first glance, figure 9 appears to be a shotgun pattern of data 
points, with transition Reynolds numbers varying from about 1 million to 30 million. However, 
even within this jumble of data there are some definite trends. First, note that the flight data give 
the highest transition Reynolds numbers and the wind tunnel data give the lowest values. The 
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ballistic range points fall, more or less, in the middle. The tendency to predict early transition in 
many wind tunnels has been widely observed 36 - 37 and correlated by Dougherty and Fisher 37 with 
the intensity and, to some extent, the frequency of the disturbances in the facilities. The problems 
encountered in wind tunnel test measurements of transition have encouraged the use of ballistic 
ranges. 38 - 39 More recently, the development and use by Beckwith 40 of a low-disturbance, high 
Reynolds number, supersonic wind tunnel has yielded valuable data. 41 - 42 For example, Chen 
et al. 42 measured similar transition Reynolds numbers at a Mach number of 3.5, on flat plates and 
slender cones with adiabatic walls. In high-noise wind tunnel tests, the transition Reynolds num- 
bers on flat plates were only half as large as the values measured on cones, 42 implying that the flat- 
plate boundary layer was more sensitive to free-stream disturbances. 

Empirical Correlations 

Some correlations for the beginning of transition found in the open literature will now be 
discussed. Despite the data scatter in figure 9, there is a discernible trend of increasing transition 
Reynolds number with rising Mach number. The same trend was observed by Softley et al. 43 and 
is shown in figure 10 (taken from Ref. 43).The data in figure 10 are for sharp cones, having 
Reynolds numbers based on nose radii of less than 200. The high stability of the laminar boundary 
layer to disturbances at hypersonic edge Mach numbers has been observed by other researchers 44 
and lends credence to the trend shown in figure 10. However, the influence of wall temperature on 
hypersonic transition appears to be an open issue, with some tests showing strong effects and 
other tests indicating little effect. 43 

Unlike the hypersonic boundary layer, the subsonic one existing on the blunt noses of 
high-speed flight vehicles is easily tripped by surface roughness. The mechanism has been exten- 
sively studied and correlations have been published 45-47 Since strong pressure gradients exist on 
the blunt noses, the correlations are for Reynolds numbers based on boundary layer momentum 
thickness rather than the body lengths used for sharp cones and flat plates. A formula suggested by 
Laderman 46 is 

m fir T 

Re«bt = 215 / (g-T^J (50) 


where k/0 is the ratio of roughness element height to local momentum thickness. (Charts of 
momentum thickness in high-speed flight for bodies with various amounts of nose bluntness can 
be found in Ref. 48.) Another correlation for surface-roughness-induced transition is presented by 
Amirkabirian et al 49 for the Shuttle orbiter and is shown in figure 11. Although the Shuttle tiles 
are very smooth, the "roughness" results from misaligned tiles and the gaps between tiles. Again, 
most of the flight data points are well above the shaded band, which is based on wind-tunnel tests. 
An approximate correlation is 49 


Ree/Me = const. (51) 

where the constant varies from 150 to 350, depending on the ratio of roughness height to momen- 
tum thickness, etc. Additional transition data are available for configurations of practical interest 
such as blunted cones at angle of attack 50 and on swept- wing leading edges. 51 Despite the high- 
transition Reynolds numbers which can be expected at hypersonic speeds, full-sized vehicles will 
still experience turbulent boundary layers over much of their surface. The heat transfer resulting 
from turbulent boundary layers will be discussed next. 



Turbulent Heat Transfer— Cold Wall 
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For the laminar boundary layer, it was shown that exact solutions could be found which 
model the viscous flow over important regions of a vehicle. In contrast with the laminar boundary 
layer, there is no complete theory for modeling the turbulent transport of mass, momentum, and 
energy. The equations for the turbulent boundary layer can be formulated only with the aid of 
experimental data such as mixing lengths, etc. Therefore, all turbulent "theories" are semiempirical 
and care must be taken when these formulations are applied beyond the range of conditions for 
which they were derived. 

Comparison of Approximate Methods 

A number of approximate methods have been widely used to calculate the turbulent bound- 
ary layer heating. All of the methods rely on Reynolds analogy which relates skin friction and heat 
transfer. 8 In addition, all methods use some form of coordinate transformation for extending the 
well-verified incompressible skin-friction and heat-transfer formulas to the compressible, high- 
speed conditions. The transformations therefore are functions of Mach number, Reynolds num- 
bers, and temperatures at the boundary layer edge and at the wall. Among the most widely used 
approaches are the reference enthalpy, or reference temperature, Spalding and Chi, Van Driest n, 
and Coles. These methods have been extensively compared and their formulations explained in the 
literature by Hopkins and Inouye 52 and by Cary and Bertram 53 ; therefore, only brief descriptions 
will be presented here. 

The reference temperature concept was originated by Rubesin and Johnson 54 and the coef- 
ficients modified by Sommer and Short. 55 The basic assumption is that a temperature within the 
boundary layer can be calculated, using an expression with empirically determined coefficients, 
which will yield the correct skin friction in a compressible, supersonic boundary layer when the 
incompressible relations are used. Although the method was first formulated for laminar boundary 
layers, it was found to give good results for turbulent boundary layers also. Eckert 56 extended the 
method to high-speed flows in equilibrium having real gas effects by using enthalpy in place of 
temperature. Van Driest II refers to Van Driest’s second method which uses the von Karman 
mixing length. 57 In addition, the Crocco temperature profile is assumed through the boundary 
layer. The Spalding and Chi method 58 uses Van Driest's formulation to relate the incompressible 
and compressible boundary layer skin friction. Empirical expressions are used containing the 
boundary layer edge, wall, and recovery temperatures. The correlation functions were empirically 
determined using a large body of data. However, some of the data contained systematic errors yet 
all data were weighed equally. The method of Coles 59 is based on the concept of a constant mean 
temperature of a layer, called the substructure, which is located near the wall but extends beyond 
the laminar sublayer. The substructure temperature is based on the boundary layer edge Mach 
number and temperature and the wall temperature and is used to evaluate a Reynolds number. 
Additional mathematical details on the above four methods can also be found in Ref. 8. 

The ability of the four methods to predict turbulent heat transfer has been compared by 
several authors. The comparisons of Hopkins and Inouye 52 are shown in figure 12 for the Mach 
number range 4.9 to 7.4 and for ratio of wall temperature to recovery temperature from 0 1 to 
0.82. The conclusion drawn from figure 12 and comparisons in Ref. 52 is that the Van Driest II 
method gave the best overall results while Spalding-Chi consistently underpredicted the heating. 
However, note that disagreement between the prediction and the data of 20% to 30% is common. 
The conclusions of Ref. 52 are corroborated by the results of Chien 60 in figure 13 for a wall tem- 
perature ratio of 0.2 or greater. Chien’s tests were conducted at Mach 7.9. For a very cold wall 
temperature ratio of 0.1 1, the Spalding-Chi method is shown (fig. 13) to be best. However, shock 
tube experiments 61 - 62 at wall temperature ratios from 0.01 to 0.24 and at Mach numbers of 1.5 to 
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3.1 indicate the best agreement with Van Driest II in Ref. 61, while Ref. 62 supports using 
Spalding-Chi. In fact, Ref. 53 and Cary 63 concluded that the Spalding-Chi method yielded the 
best results over the entire temperature ratio range 0. 1 to 0.7 at Mach numbers of 4 to 13. Charts of 
turbulent heat-transfer coefficients based on the Spalding-Chi correlations were made by Neal and 
Bertram. 64 Zoby and Graves 65 statistically studied the effect of using two different forms of 
Reynolds analogy on the Van Driest II, Spalding-Chi, and reference enthalpy methods. Data from 
ground-based facilities and from flight tests were analyzed. It was concluded the Colburn modifi- 
cation of Reynolds analogy, which is 


St = 


Cf 

2Pt0.667 


(52) 


usually resulted in somewhat better agreement with test data than using von Karman's version of 
Reynolds analogy. A comparison of the three heat-transfer prediction methods with the data 65 
typically yielded rms errors of 15% to 18%. 

The conclusion is that all the semiempirical methods give only approximate results that can 
be off by 20%, or sometimes more. However, for wall temperature ratios below 0.2, Spalding-Chi 
does seem to predict the heating better than the other three methods. It should also be noted that 
finite-difference codes for calculating perfect gas turbulent boundary layer properties have been 
developed by Wilcox 66 " 68 and, no doubt, others. 

Reference Enthalpy Method 

The reference enthalpy method has been widely used for three decades. The most common 
expression for the reference value of the enthalpy 56 is 

h' = 0.22h aw + 0.28h e + 0.5h w (53a) 

although a number of variations have been used. At least one other version 8 given by 

h' = 0.18h aw + 0.32h e + 0.5h w (53b) 

yields slightly better results. The Eckert version, Eq. (53a), was used by Arthur et al. 69 to derive 
simple, closed-form correlation equations for the turbulent flat-plate heating. The Blasius incom- 
pressible, turbulent skin-friction relation was used in Ref. 69 for Reynolds numbers to 10 million 
in the form 


cf 0.0296 
2 (Re') 0 - 2 


and the Schultz-Grunow equation 


c f 0.185 
2 (logio Re') 2 - 584 


(54) 


(55) 


for Reynolds numbers above 10 million. In Eqs. (54) and (55), the prime denotes that the quan- 
tity is evaluated at the conditions corresponding to the reference enthalpy. The thermodynamic and 
transport properties of Refs. 5 and 6 were used to calculate the heat-transfer rate at altitudes from 
3 to 60 km and at velocities from about 1 to 1 1 km/sec. The resulting correlation, which is valid 
to ±15% for Voo > 1500 m/sec, but less than 3960 m/sec 
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In Eq. (56) the Newtonian approximation for pressure and velocity were used in the form 

Pe “ P~V* sin 2 6 (57) 


u e = Voo cos 5 (58) 

Therefore, the heat transfer is underpredicted at small values of 5. (For example, at 5 = 10°, the 
heating will be about 6.5% too low, which is wellwithin the uncertainty of the reference enthalpy 
method.) Note that for the beginning of the turbulent boundary layer, the location at which transi- 
tion begins, x bt , was chosen. While this choice is somewhat conservative, it is supported by the 
correlations of Refs. 53 and 65. The turbulent heat transfer on a sharp cone of half angle, 5, 
would be given, approximately, by 8 

qw c =1.15q WF p (59) 

The relations derived so far are valid if the boundary layer is either laminar or turbulent. In 
locating the effective origin of the turbulent boundary layer in Eq. (56), it was assumed that the 
transitional flow region is equal in length to the preceding laminar flow distance. For lack of a bet- 
ter model, it is assumed here that the heating rate in the transitional region varies linearly with dis- 
tance, having the laminar value at the beginning of transition and the fully turbulent value at the end 
of transition. This approximation gives, roughly, the correct average skin friction according to 
Dhawan and Narasimha, 70 and, therefore, the correct heating rate. 

It is now possible to calculate the heat-transfer rate on a simple configuration. The example 
chosen here is the heating along the fuselage windward centerline of the Shuttle orbiter during the 
first flight (see fig. 14) from Tauber and Adelman. 71 The calculations were made using Eqs. (48) 
and (56) by assuming that the bottom surface was a wedge at the local surface slope and in 
radiative equilibrium. The agreement between the measured values and the calculations is fairly 
good. The turbulent calculations are definitely too high, possibly because of the cold wall condi- 
tions or the conservative choice of the origin. The error bars on the data represent the effect of 
uncertainties in the surface emissivity. 72 At the flight conditions shown in figure 14, molecular 
dissociation was sufficiently small to make catalytic wall effects unimportant in the laminar bound- 
ary layer. Finite-rate wall catalysis has not been observed in turbulent boundary layers, probably 
because of the intense mixing. 

Cross-Flow and Surface Roughness Effects 

Three-dimensional shapes, or axisymmetric bodies at angle of attack, all experience cross 
flows. The heating of such complex flows must, generally, be computed using finite-difference 
methods. 73 ’ 74 However, for a few geometrically simple shapes, or in the plane of symmetry, 
approximate methods can be used. One such example is the calculation of the turbulent heating on 
the leading edge of a swept wing. The cross flow induced by the leading-edge sweep promotes 
boundary layer transition. 31 Also, a turbulent boundary layer at the wing-fuselage junction can 
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cause transition at the leading edge. To approximately account for the turbulent flow along the 
leading edge, Eq. (56) must be used for the flat-plate contribution in Eq. (49). 

Another example where the heating in a complex flow can be treated approximately is the 
sharp cone at angle of attack. For this case, the tangent cone approximation® yields reasonably 
good results when the angle of attack is less than the cone half angle. In the tangent-cone method, 
an imaginary cone is used having a half angle which equals the slope of the local ray on the cone at 
angle of attack. However, when the cone also has a blunt nose, the heating distribution becomes 
more complex (see fig. 15, from Widhopf 75 ) and must be calculated numerically. 74 

Vehicles designed to fly at hypervelocities are usually protected by ablative heat shields. 
Frequently, ablation causes a rough surface which can trigger early boundary layer transition. 
Roughness-induced transition on blunt noses of bodies, where the boundary layer is thin and the 
flow subsonic, has been widely observed and studied. 45-47 In addition to triggering premature 
transition, the surface roughness can increase the turbulent heating by over 100%. The turbulent 
heating caused by roughness elements equal in size to about 1% of the nose radius is illustrated in 
figure 16 (from Chen 76 ). The peak rough- wall turbulent heating formulation of Chen predicts 
more than twice the smooth- wall value and agrees reasonably well with the heavily scattered data. 
Note, also, that the peak heating rate occurs about 30° off the stagnation point, compared to 35° for 
the smooth surface value. The computation is performed by numerically integrating a set of modi- 
fied turbulent boundary layer equations which are based on the assumption that the effect of 
roughness on the mixing is restricted to a relatively thin layer near the wall. 

Next, the heating associated with regions of separated flow will be discussed. This will be 
followed by a brief coverage of shock-interaction heating. 

Separated Flow and Shock-Interaction Heating 

Although regions of separated flow occur frequently when control surfaces are deflected to 
large angles on hypersonic vehicles, the computations of such viscous flow fields is complex and 
relies heavily on definite-difference, numerical methods. Another mechanism that can cause flow 
separation, and severe local increases in heating near the flow reattachment point, is the presence of 
gaps, or cavities, and steps in the surface. The discussion of heating in cavities will be followed by 
a brief coverage of the heating experienced in the base regions of bodies. The extremely high local 
heating that can occur when shocks intersect is covered last. 

Shock-Induced Separated-Flow Heating 

In a supersonic or hypersonic flow, a shock is generated when the local surface slope 
increases abruptly. For example, such a change in slope occurs when a control surface is deflected. 
Since the static pressure increases through the shock, the higher pressure is fed upstream through 
the subsonic part of the boundary layer, which causes the flow to separate (see fig. 17). For a 
given pressure increase, the extent of the separation region in a laminar boundary layer is greater 
than in a turbulent one since the former has a greater fraction of subsonic flow. A region of recir- 
culating flow is formed in the comer with the fluid adjacent to the wall moving upstream. The 
heating in the separated-flow region is increased over the value ahead of the separated region and 
peaks at, or slightly beyond, where the boundary layer reattaches on the compression surface. The 
heating in separated and reattached flows has recently been reviewed by Merzkirch et al. 77 

Extensive heat-transfer measurements were performed by Holloway et al. 78 in laminar, 
transitional, and turbulent separated flows. The experiments were performed at Mach 6 and at 
maximum length Reynolds numbers from 0.9 million to 7.3 million. An example of the heat 
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transfer in a turbulent boundary layer at Mach 6 with attached and separated flow on a 75° swept- 
back, sharp delta wing is shown in figure 18 (from Keyes et al. 79 ). The wing is a t zero angle of 
attack. The flap deflection of 20° does not cause flow separation, whereas the deflection of 40° 
produces strong separation. (The problem of trying to calculate the heating x on the flap is evident 
from the variety of methods used.) Various researchers have correlated the ratio of turbulent heat- 
ing rates ahead of the separated region to the peak rates on the flap with the corresponding pressure 
ratio. Back and Cuffel 8 ® give the empirical relation 

ch 2 /chi = (P2/Pi) 085 (60) 

where their (dimensional) heat transfer coefficient is defined as 


Ch = q/(h aw - hw) (61) 

Nestler 81 suggests an exponent of 0.80 in Eq. (60) and Alzner and Zakay 82 find values from 
0.725 to 0.815 for turbulent flows. 

For hypersonic laminar boundary layers, Mathews and Ginoux 83 suggest a relation in 
terms of Stanton number at peak heating, St m ax. which is approximately 

Stmax ~ (Rei)-°- 5 (62) 

where Rei is the Reynolds number ahead of the separated region on the flat plate. The constant of 
proportionality appears to be a strong function of wall temperature ratio. 

The problem of determining the heating when flow separation occurs has been attacked 
successfully numerically. Shang and Hankey 84 calculated the supersonic flow field for the turbu- 
lent boundary layer case, while Hung and MacCormack 85 did the hypersonic, laminar case. More 
recently, Lawrence et al. 86 used a perfect gas parabolized Navier-Stokes (PNS) code to calculate 
the hypersonic flow on a sharp flat plate at zero angle of attack with a flap (fig. 19). (The flow 
field shown in fig. 19 could simulate the hypersonic flow over a sharp-nosed airfoil with a large 
flap.) A comparison of calculated laminar heat transfer with experimental data from Ref. 86 is 
shown in figure 20. At a Mach number of 14, the 15° deflection angle did not cause significant 
separation. A newer version of the code, incorporating real gas effects, was used to perform the 
same calculation as above, but at twice the previous Reynolds number and at a Mach number of 
20- 7 The result is shown in figure 21. (Since the calculations were performed on different com- 
puters, no direct comparison of computation time is possible.) The heat-transfer coefficient of 
Refs. 86 and 87 is defined as Ch = Pr St(p e Ue/pooV«). 

Shock-induced boundary layer separation and increased local heating also occur when a 
shock wave impinges on the flow over a flat plate. Experiments performed by Back and Cuffel 88 
m a turbulent boundary layer at Mach 3.5 confirmed the empirical relation given by Eq. (60), but 
gave an exponent of about 0.73. Skebe et al. 89 performed detailed boundary layer measurements 
of the effects of shock interactions on a flat plate having initial laminar, transitional, and turbulent 
flow. The Mach numbers ranged from 2 to 4 and Reynolds numbers per meter varied from about 
4.7 million to 30 million. Detailed boundary layer surveys were made, but heat transfer was not 
measured. Hodge 9 ® compared numerical calculations using MacCormack's code 85 with pressures 
measured for laminar boundary layer shock interactions at Mach 8. He found good agreement 
except near the beginning of the induced separation region. Hung et al. 91 present correlations for 
the peak heating from shock impingement on a flat plate at hypersonic speeds. The following 
expressions for the Stanton numbers at peak heating are given in Ref. 91: 
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for laminar/laminar flow 

St m ax - 0.664(Re , x )-^ 5(Pr , )-0-667 (63) 

for laminar/turbulent and turbulent/turbulent flow 

Stmax - 0.0576(Re , x )“°- 2 (Pr') -0 - 667 (64) 

where the prime denotes that the parameter is based on the reference enthalpy value (Eq. 53a) at 
the peak heating location. Viegas and Coakley 92 studied turbulence models that were used in the 
Navier-Stokes equations to calculate shock- separated turbulent boundary layers. They concluded 
that the one-equation model proposed by Glushko and Rubesin gave superior results, including 
more accurate skin-friction values, than the zero equation model. 

Heating in Cavities 

Cavities occur on all vehicles and can cause large increases in local heating at high speeds. 
First, cavities can trigger boundary layer transition. Second, the heating within the cavity can be 
much higher than the heating of the surface just ahead of the cavity. Experiments of heating within 
cavities in turbulent flows were reported by Nestler et al. 93 and some of their results are shown in 
figure 22. The ratio of the flat-plate heating rate to that measured within the cavity is shown in fig- 
ure 22; Ch is defined by Eq. (61). (The curves have been drawn through the original data points, 
which are deleted.) The Mach number on the plate was 6.3 and the local Reynolds number was 
about 6 million. The turbulent boundary layer thickness just ahead of the cavity location was 
1 cm. For the cavity length-to-depth ratio (L/H) used in figure 22 of 5, 15, 30, the cavity depth 
was 2.5 cm, 2 cm, and 1 cm, respectively. Note that for L/H = 5, the flow passes over the 
cavity and the heating on the cavity floor is decreased from the undisturbed flat-plate value. When 
L/H » 1 1, the nature of the flow in the cavity changes according to Gortyshov et al. 94 For values 
of L/H larger than 1 1, the flow expands into the cavity, attaches on the floor, and separates near 
the downstream comer, where the heating near reattachment reaches three times the undisturbed 
value. 


For supersonic turbulent flow, the average heat transfer in the cavity is 94 

Stave = 0.48(Re L )-°- 4 (L/H)0.2 (65) 

where the Stanton and Reynolds number properties are evaluated using the undisturbed flow con- 
ditions and the cavity length. 77 For cavities in supersonic and hypersonic laminar boundary layers, 
Lamb 95 gives the following expression: 


St - 2.07 (Res) -1 (66) 

In Eq. (66), the Reynolds number is, again, evaluated at the undisturbed flow conditions and is 
based on the boundary layer thickness ahead of the cavity, while the Stanton number is related to 
the heating on the leeward surface. 77 

Heating in Base Regions 

Since the heating in the base regions of hypersonic vehicles is one to two orders of magni- 
tude less than on the windward portions of most vehicles, the topic will be discussed only briefly. 
Francis 96 presents flight data for ablating very slender (4.5° half-angle) cones with turbulent 
boundary layers. For a blunted cone, the heating in the base region was, approximately, 6 % of the 
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cold wall value without ablation on the cone afterbody. Bulmer 97 - 98 gives a correlation for turbu- 
lent flow base heating on cones in hypersonic flow in the form 

35.5(Ret>/Re c ) 2 - 2 (67) 

In Eq. (67), the subscript b refers to the evaluation at conditions where the flow has expanded 
isentropically to the base pressure and the reference length is the base radius. The subscript c 
refers to the flow conditions on the cone just ahead of the base. 77 For base heating where a laminar 
boundary layer exists. Ref. 77 suggests using Eq. (66). 

Shock-Interference Heating 

Shock-interference flows occur when, for example, the bow shock from the blunt nose of a 
vehicle impinges on the leading edges of wings, fins, and inlet cowl lips. The analysis of shock- 
interference heating is fundamentally complex, but cannot be ignored because heating rate increases 
of one order of magnitude have been observed. The intersection of strong shock waves can pro- 
duce free shear layers, jets, or expansions. Different shock interactions and the resultant flows 
- were delineated and studied by Edney." In Ref. 99 it was found that large increases in heating 
result from free shear layer attachment and supersonic jet impingement. Keyes and Morris 100 per- 
formed experiments at a free-stream Mach number of 6. Turbulent shear layer attachment on a 
hemisphere was found to produce local heating rates which were up to 14 times higher than the 
undisturbed stagnation point value. An approximate correlation for the ratio of peak local heat 
transfer to the (undisturbed) stagnation point value 100 can, again, be written using the corre- 
sponding pressure ratio, in the form 


Chp/Ch s - C(p p /p s ) 1 -35 (68) 

where C = 1.58 when the shear layer is turbulent and C = 0.63 when it is laminar. The 
dependence of the heating on the shear layer state (laminar, turbulent, or transitional) is, again, 
emphasized by Keyes. 101 It is shown that the initial state of the shear layer can have a major influ- 
ence on the heating. 


Heat-Protection Methods 

The preceding discussion has been limited to methods of calculating the heating rates at the 
body's wall in the absence of interactions with the surface, i.e., the cold wall heating. The follow- 
ing discussion will deal with methods of protecting the body surface from the aerodynamic 
heating. 

Several heat-protection methods are widely used. If the heating rate and duration of the 
pulse is modest, a high-heat capacity metal heat sink, such as Beryllium oxide, can absorb about 
6.3 MJ/kg (about 2700 Btu/lb). For many purposes, a more effective mechanism of heat absorp- 
tion is to transfer mass from the surface to the boundary layer. The three main methods of mass- 
addition cooling are described below. 

Film cooling: fluid is injected into the boundary layer near the stagnation point, from 
where it spreads back over the body. 

Transpiration cooling: the fluid passes through a porous wall into the boundary layer. This 
differs from film cooling in that fluid injection can occur all along the wall. 
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Ablation: the surface material is allowed to melt and vaporize, and enters the boundary 
layer in a liquid or gaseous state. The mass addition of fluid into the boundary layer reduces heat 
transfer by (1) absorbing heat through phase changes of the material from solid to liquid to gas 
and (2) thickening the boundary layer and altering the velocity profile. Ablative heat shields are 
frequently impregnated with carbon fibers. The fibers provide structural reinforcement and, at high 
temperatures, form a char layer. The char is porous and permits percolation of gases and reradiates 
heat. For example, at a temperature of 3000 K, over 400 W/cm can be reradiated. 

To protect the structure and the interior components of the vehicle, insulation material is 
required between the heat shield and the structure. Since the heat transfer occurs by conduction, 
materials having low conductivity and low density are used. 

Mass- Addition Cooling 

The most widely used method of protecting vehicles against the high heating encountered in 
hypervelocity flight is mass-addition cooling. An early study of (stagnation region) film cooling 
was made by Howe and Mersman. 102 Low 103 calculated the effect of transpiration cooling on the 
laminar boundary layer on a flat plate in supersonic flow where the injectant was air. Subse- 
quently, Pappas and Okuno 104 ' 1 ^ measured the influence of injecting various gases, which were 
both lighter and heavier than air, on the heat transfer on cones, having laminar and turbulent 
boundary layers, in supersonic flow. Kaattari 107 and others performed mass-injection experiments 
on blunt-nosed bodies, including hemispheres. The results of the fluid-injection experiments were 
correlated for laminar boundary layers in Ref. 23. For the turbulent case, the coefficients were 
found empirically. The correlation is, using an expression of the form 7 

\\f = q w /qwB=o = 1 “ a(me/fn w ) n B + b(ffvfh w ) 2n B 2 (69) 

where rhe and m w are the molecular weights of the shock-layer gas and the injection gas, respec- 
tively. The dimensionless mass-addition parameter, B, is defined as 


B= J^(He-hw) 

9wb=0 


(70) 


where m w is the mass-addition rate. The constants a, b, and n in Eq. (69) have the following 
values: 


a b n 


Stagnation point (lam.) 0.72 0.13 0.25 

Cone, or flat plate (lam.) 0.795 0.11 0.25 

Cone, or flat plate (turb.) 0.29 0.015 0.5 

The results from Eq. (69) are compared with the measurements of Refs. 104-107 by plot- 
ting the heating reduction as a function of the mass-addition parameter. In figure 23, the heating 
reduction is shown when air is injected into an air boundary layer, m w = me, at the stagnation point 
and on cones, or flat plates, having laminar and turbulent flow. The agreement is good. It is also 
apparent that a much higher injection rate is required to cool the wall when turbulent flow exists, as 
confirmed by Dershin et al. 1 ^ The effect of molecular weight ratio on the agreement between 
Eq. (69) and the data is illustrated in figure 24 for laminar flow. Again, the agreement is reason- 
ably good. The beneficial effect of injecting gases which are lighter than air is evident. Since lighter 
gases thicken the boundary more than heavy ones, the velocity gradient in the boundary layer is 
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reduced more. (For Freon, m w ~ 120 was used.) The influence of molecular weight ratio is more 
pronounced for the turbulent boundary layer than for the laminar one, as is evident in figure 25. 
For gases that are heavier than air, Eq. (69) underpredicts the wall cooling. Note also in figure 25 
that in the turbulent flow case, the mass-injection efficiency decreases with increasing edge Mach 
number. This observation is confirmed by Jeromin 109 in his lengthy survey paper and is shown 
more dramatically in figure 26 (based on Ref. 109). 

Expressions such as Eq. (69) can be used to approximate the behavior of relatively low- 
temperature vaporizing ablators, for example, Teflon. For high-temperature, carbonaceous abla- 
tors, the presence of chemical reactions between the injected and resident species in the boundary 
layer must be included in addition to the surface radiation. Neither of these phenomena is 
accounted for in Eq. (69). The performance of carbon ablators was calculated by Putz and 
Bartlett, 110 who derived empirical relations, similar to Eq. (69), for several materials. A stagna- 
tion point correlation for graphite, from numerical solutions, is 110 

\jr' = 1 - aB + bB 2 + cB 3 - dB 4 (71) 

where a = 0.6563, b = 0.01794, c = 0.06365, and d = 0.01 125. The previously defined 
heat- transfer ratio, \\r, is related to" Eq. (71) by 

y=v ' + ^)/ ( 1_h *) ™ 


where the enthalpy ratio AHc/He accounts for gas-phase chemical reactions and diffusion effects 
for injected species. The variation of \| f with the mass-addition parameter, Eq. (71), is shown in 
figure 27 and the enthalpy ratio from chemical reactions is given in figure 28. (Both figures are 
from Ref. 110.) 

Slender vehicles can experience very high heating at the nose in hypervelocity flight, caus- 
ing surface recession and shape change at the nose. If the nose ablation is asymmetric, an undesir- 
able pitching moment results. The slender-vehicle nose-ablation problem is formulated and the 
results of numerical solutions are presented by Chin. 111 The roughness of the ablating surface can 
also cause premature boundary layer transition, which further intensifies the heating, as previously 
mentioned. The effects of nose surface roughness on transition and heating are discussed by 
Grabow and White. 112 Slender-body nose-shape changes as determined from wind tunnel tests 
using low-temperature ablators were measured by Kobayashi and Saperstein. 113 Using film cool- 
ing to reduce the surface recession and shape change of slender noses was studied by Gold 
et al. 114 The heat protection by mass injection was also theoretically studied for other geometric 
shapes and will be briefly covered next. 

The effect of mass injection on the laminar heating along the stagnation line of an infinite, 
swept cylinder was calculated by Libby and Kassoy. 115 For massive blowing, it was found that 
the spanwise flow separated, but the chordwise flow remained attached. Wortman analyzed the 
effect of mass addition on several blunt shapes ranging from hemispheres to cylinders. 11 6,H7 j n 
Refs. 116 and 1 17 the laminar boundary layer equations are solved for a real gas in equilibrium 
and for a variety of light and heavy (compared to air) injected gases. The results are presented in 
the form of correlation functions. In Ref. 118, Wortman calculates the reduction of laminar heat- 
ing by mass injection on sharp cones at angle of attack. It was found that normalizing the heat 
transfer by the zero angle of attack values eliminated the influence of injected gas properties. Inger 
and Gaitatzes 1 19 computed the effect of strong blowing on high Reynolds number laminar flow, at 
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supersonic speeds, about slender bodies, including induced pressure effects. Some of the effects 
computed in Refs. 118 and 119 are open to question, however. For example, it is well known 
that cross flow on cones at large angles of attack 118 can induce boundary layer transition. It has 
also been observed that massive blowing at high Reynolds numbers can cause transition very close 
to the stagnation point of a blunt body, as will be discussed next. 

Iniection-Induced Transition 

Demetriades et al. 1211 and Kaattari 1 ^ observed boundary layer transition caused by mass 
injection on blunt bodies in hypersonic flows. The effect was observed at Reynolds numbers based 
on body diameter as low as about 1 million. Similar results were observed by Wimberly et al. 121 
on a slender cone tested at a free-stream Mach number of 12.2 and a free-stream Reynolds number 
of 4 million. In most of the cases observed, the occurrence of transition raised the heating above 
the laminar value without mass addition. Park 122 presented calculations to explain the injection- 
induced early transition on blunt bodies. However, no general method is available for predicting 
the injection rates that will lead to premature transition. 

Ablation Material Response 

To calculate the response of an ablation material to high heating rates, it is necessary to 
know the high-temperature properties of the material. Eighty-five different graphitic materials were 
tested by Lundell and Dickey. 12 3,124 The pressures ranged up to 4.5 atm and temperatures up to 
4000 K. The materials that were studied included ATJ graphite, both two- and three-dimensional 
carbon-carbon composites, pyrolytic graphite, mesophase graphite, glassy carbon, and natural 
graphite. It was found that ATJ graphite performed slightly better than the single-phase materials 
(pyrolytic graphite, mesophase graphite, and glassy carbon). However, at temperatures above 
3600 K the ATJ mass loss increased exponentially with temperature. The mechanism responsible 
for the greatly increased mass loss was particulate removal, or spallation, probably caused by 
compressive thermal stress. 

A different class of ablation materials consisting of ceramics was studied by Ziering. 126 
Metallic and ceramic materials have much greater mechanical strength than graphite, for example, 
and would experience much less erosion when traversing dust or rain clouds. Tests of the follow- 
ing materials were reported in Ref. 125: silicon nitride, silicon carbide, tantalum carbide, and 
tungsten. Silicon nitride's mass loss, surface recession, and energy absorbed per unit mass were 
found to be by far the best of the group, with silicon carbide rating a distant second. The energy 
absorbed per unit mass ablated was 6 MJ/kg (2600 Btu/lb) for silicon nitride, while silicon car- 
bide's value was 0.72 MJ/kg (310 Btu/lb). 

For additional information on the application of mass-transfer cooling, Hartnett 126 is 
recommended. 


Einite-Difference Codes 

A large number of computer codes exist which can be used to calculate heat transfer at high 
velocities. Most of these codes were not intended for widespread use and have never been docu- 
mented. By necessity the following discussion will be limited to a small, but representative, sam- 
pling of codes that are being used to calculate cold- wall aerodynamic heating. The discussion con- 
cludes with a description of some codes that can be used to compute surface ablation. 

Thompson et al. 122 compared the heating on blunted slender cones at angles of attack as 
calculated by four different codes with both wind tunnel and flight data. The codes consisted of a 
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three-dimensional, viscous flow, finite-difference method called VSL3D, and three approximate, 
engineering- type codes known as MINIVER, AEROHEAT, and INCHES. As expected, VSL3D 
gave the best agreement with the data. Of the approximate methods, MINIVER predicted the nose 
bluntness effects poorly and AEROHEAT underpredicted the laminar flow heating at angle of 
attack by as much as 30%. The INCHES code gave good results at zero angle of attack, but at 
small angles of attack the windward ray heating was underpredicted by up to 40%. However, the 
engineering-type codes all used small amounts of computer time and their utility must be evaluated 
in that light. 

Space-Marching Codes — PNS 

The most widely used methods to calculate hypersonic flow fields and heat transfer are the 
downstream marching solutions of the parabolized Navier-Stokes equations known as PNS codes. 
The parabolized approximation is valid if the steady-state shock-layer flow, which must be super- 
sonic, and the subsonic portion of the viscous-layer flow are both in the streamwise direction. 
Therefore, flows with large amounts of separation, and the attendant reverse flow, cannot be 
treated; however, cross-flow separation is permitted. When the above limitations are imposed, the 
Navier-Stokes equations becomes parabolic in the streamwise direction, thus permitting the down- 
stream marching from an initial, "supersonic data plane. Chaussee 128 evaluated the ability of a per- 
fect gas PNS code to predict heat transfer on a blunted bicone having laminar and turbulent bound- 
ary layers. Comparisons with Mach 8 wind tunnel data gave reasonably good results when the 
damping coefficients (artificial viscosity) were held at the smallest value giving stable solutions. 
Chaussee and Rizk 129 used the same PNS code to make calculations of the flow over deflected 
control surfaces on several body shapes. The flap deflections were kept small enough to avoid 
separating the turbulent boundary layer. Although heating rates were not calculated, pressure dis- 
tributions and other flow-field information were presented. Rizk et al. 130 used the PNS code to 
calculate heating on a bicone at angle of attack on the windward, leeward, and 90° plane. Their 
results were compared with wind tunnel measurements made at Mach 10 and are shown in fig- 
ure 29 (from Ref. 130). The comparison shows some disagreement near the nose of the vehicle. 
The differences on the aft part of the lee side are due to boundary layer transition which the code 
cannot predict. 

The PNS code calculations of Refs 128-130 were made assuming a perfect gas. It would 
be surprising if the aerodynamic heating in a Mach 10 flow could be accurately predicted using the 
perfect gas assumption. In fact, Balakrishnan and Chaussee 131 showed that using perfect gas cal- 
culations caused the stagnation point heating to be significantly underpredicted at high speeds. 
Efforts to explain Shuttle orbiter surface-temperature measurements provided one impetus to 
include real gas effects in PNS codes. Prabhu performed this task and added the capability to com- 
pute chemical nonequilibrium effects. 132 - 133 Subsequently, Prabhu et al. 134 extended the PNS 
code to permit three-dimensional flow computations. Examples of Prabhu's calculations are shown 
in figure 30, for a real gas in chemical equilibrium, and in figure 31 for the chemically reacting 
flow. Both figures 30 and 31 are for sharp cones having a 10° half angle, flying at a speed of 
8100 m/sec at an altitude of 61 km. The small cone angle results in an equilibrium shock layer 
temperature of 3100 K at zero angle of attack. The relatively low temperature limits the 
nonequilibrium effects. Tannehill et al. 135 modified a PNS code 86 to include chemical reactions 
and also added an upwind-differencing algorithm. The upwind-differencing method improves 
shock capturing and avoids the instabilities inherent in central-differencing schemes caused by 
flow-field discontinuities such as shocks. The amount of artificial viscosity that is needed to con- 
trol the oscillations occurring in central-differencing schemes can lead to inaccuracies in the com- 
puted results, as was observed in Ref. 128. 
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Time-Iterative Codes 


To overcome the inherent limitations of the PNS codes, new methods are being developed 
to calculate the hypervelocity, real gas flow about three-dimensional bodies. The new codes use a 
time-iterative procedure which eliminates the need for artificial viscosity. However, the computa- 
tion time for the time-dependent methods is nearly two orders of magnitude greater than for the 
downstream marching codes. 

One example of a recently developed time-iterative code is the work of Palmer. 1 36 Palmer's 
code employs an implicit, flux- splitting, shock-capturing method and includes equilibrium, real- 
gas effects. The code is being extended to treat chemical reactions. 137 

Edwards et al. 138 compared four computer codes by using each to calculate hypersonic 
flow over two different bodies. The four codes consisted of an upwind PNS solver, UPS, 86 and 
three time-iterative methods which solved the three-dimensional, thin-layer Navier-Stokes equa- 
tions. The time-iterative codes consisted of a central-differencing, time-variational, diminishing 
code known as VAIR3D/TVD, a partial flux-split code called F3D, and an upwind solver, UWIN. 

A comparison of wind tunnel measured heating rates with calculations using the above codes is 
shown in figure 32 (from Ref. 138). The comparison is for a bicone model at Mach 10 and a 10° 
angle of attack. All the codes give similar results with minor exceptions. Near the body's nose, the 
heating in the 90° meridional plane is underpredicted. Also, the occurrence of boundary layer tran- 
sition on the leeward side increases the measured heating well above the calculated laminar flow 
values. (None of the codes can predict transition.) A second comparison of the codes with each 
other is illustrated in figure 33. A generic, hypersonic body shape was used; the body is to be 
tested in the NASA Ames 3.5-ft hypersonic wind tunnel. Note that the heat-transfer coefficients, ~ 
computed at Mach 7.4 and a = 10°, vary by over a factor of 2. Apparently the strong cross flow 
occurring on the body at angle of attack poses difficulties for all four codes, although about 70,000 
grid points were used. Unfortunately, no experimental data are available yet to determine which of 
the four codes gives the best results under these conditions. The CPU time for the three time- 
iterative codes ranged from 37 to 90 min on the CRAY X-MP, compared to about 1 min for the 
UPS, space-marching code. Although great strides have been made using finite-difference methods 
to compute the (cold-wall) heat transfer on three-dimensional bodies, complex viscous flow fields 
still pose difficulties. 

Ablation Material Interaction Codes 


Compared to the calculation of cold-wall heating, the interaction between the boundary 
layer and an ablating surface can be much more complex. Several codes which can be used to cal- 
culate mass-addition effects with chemical reactions will be discussed. The codes can be used to 
calculate ablative heat-shielding requirements. 

The boundary layer integral matrix procedure (BLIMP) code 139 was developed to calculate 
the effects of chemical nonequilibrium in laminar 'and turbulent boundary layers over surfaces of 
arbitrary catalycities. The procedure is applicable to nonsimilar, multicomponent, laminar boundary 
layers with thermal diffusion and second-order, transverse curvature effects. The surface boundary 
conditions include coupling with transient ablation energy and mass balances. The BLIMP code 
has been used to analyze the surface catalytic effects experienced on the Shuttle orbiter. 1 ^ A rela- 
tively simple analysis for noncharring ablators was developed by Matting and Chapman. 140 The 
method of Ref. 140 is applicable to melting, vaporizing, or subliming surfaces, but treats only the 
stagnation point. Moss 1 1 developed a code for solving the fully viscous, ablating wall problem at 
and downstream of the stagnation point. Both convective and radiative heating can be treated. The 
coupling between the chemically reacting boundary layer and charring ablators was formulated by 
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Kendall et al. 141 The code is known as CMA. The computation is one-dimensional; although 
various body shapes can be used, the upstream history of the boundary layer is ignored. The sur- 
face ablation in gases other than air can also be analyzed; for example, rocket nozzle flows or 
flight through the atmospheres of other planets. 

In concluding this discussion of convective heating calculation methods, it should be 
emphasized that theoretical methods alone are insufficient to accurately predict the high-speed 
heating of a vehicle's surface. Codes are not yet precise enough to be relied upon solely to compute 
heating in complex viscous or three-dimensional flows. However, experiments are costly time- 
consuming, and frequently incapable of simulating flight conditions. The judicious use of verified 
codes and other theoretical procedures can greatly reduce the amount of experimental testing that is 
required. Only by using a combination of reliable theoretical methods supplemented by carefully 
conducted experiments can the heating of a vehicle be determined with a high level of confidence. 
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Appendix 

Listed below are some useful conversion factors. 


1 Btu 
1 Btu/lb 
lft 
llb 

1 slug/ft3 
1 W 

1 W/cm2 
1 W/m2 


1054.8 J 
2325.4 J/kg 
0.3048 m 
0.4536 kg 
515.2 kg/m3 
9.4805(10^) Btu/sec 
0.88076 Btu/ft 2 sec 
8807.6 Btu/ft2 sec 
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Fig. 2 Comparison of laminar convective heat-transfer parameter calculations. 
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Fig. 3 Stagnation point convective heating comparisons. 
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Fig. 4 Stagnation-point velocity gradient variation with body bluntness at Mach 4.76. 
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Fig. 5 Heat-transfer distribution on hemisphere cylinder. (From Ref. 28; reprinted with permis- 
sion of The American Institute of Aeronautics and Astronautics.) 




Fig. 6 Comparison of predicted and measured heat flux distributions on a circular cylinder normal 
to a stream. 



Fig. 7 Heat-transfer distribution on flat-nosed body. (From Ref. 28; reprinted with permission of 
The American Institute of Aeronautics and Astronautics.) 
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Fig. 8 Stagnation-point heat-transfer coefficient versus Knudsen number. V«, = 7.5 km/sec, 
r n = 2.54 cm. (From Ref. 34.) 



Fig. 9 Typical transition data — cones. 
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Fig. 10 Effect of local Mach number on transition Reynolds number — sharp slender cone — uni- 
form wall temperature. (From Ref. 43; reprinted with permission of The American Institute of 
Aeronautics and Astronautics.) 



Fig. 1 1 Flight values of the transition parameters. (From Ref. 49; reprinted with permission of 
The American Institute of Aeronautics and Astronautics.) 
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Fig. 12 Effect of wall-temperature ratio on predicted heat transfer. M e = 4.9 to 7.4; 
2Ch/Cf=1.0. (From Ref. 52; reprinted with permission of The American Institute of 
Aeronautics and Astronautics.) 
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Fig. 14 Comparison of centerline heating calculations with STS-1 flight data. (From. Ref. 71.) 
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Fig. 15. Turbulent heat-transfer distributions along the windward and leeward rays for various 
angles of attack. (From Ref. 75; reprinted with permission of The American Institute of 
Aeronautics and Astronautics.) 
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Fig. 16. Heat transfer to the hemisphere of very rough and smooth walls in a typical tunnel envi- 
ronment. (From Ref. 76; reprinted with permission of The American Institute of Aeronautics and 
Astronautics.) 



Fig. 17 Supersonic compression comer flow. 
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Fig. 18 Stanton number variation on 75° delta wing with trailing-edge flap. 
Re/m = 2.7 x 10 7 ;a = 0°. (From Ref. 79; reprinted with permission of The American Institute 
of Aeronautics and Astronautics.) 
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Fig. 20 Comparison of heat- transfer coefficients on sharp flat plate with 15° compression comer. 
(From Ref. 80; reprinted with permission of The American Institute of Aeronautics and 
Astronautics.) 
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Fig. 21 Comparison of heat- transfer coefficients on sharp flat plate with 15° compression comer. 
(From Ref. 87; reprinted with permission of The American Institute of Aeronautics and 
Astronautics.) 



Fig. 22 Relative heat-transfer coefficient measured along cavity floor. (From Ref. 93; reprinted 
with permission of The American Institute of Aeronautics and Astronautics.) 
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Fig. 23 Heat-transfer decrease at supersonic speeds for injection of air into air. (Laminar theory 
from Ref. 23; data from Refs. 104-107.) 



Fig. 24 Laminar heat-transfer decrease at supersonic speeds on sharp cone for three injected 
gases. (Theory from Ref. 23; data from Ref. 104.) 
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Fig. 27 Correlation of blowing parameter with graphite ablation boundary layer solutions. (From 
Ref. 1 10; reprinted with permission of The American Institute of Aeronautics and Astronautics.) 



Fig. 28 Variation of AHc/H e with blowing rate and total enthalpy. (From Ref. 1 10; reprinted 
with permission of The American Institute of Aeronautics and Astronautics.) 
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Fig. 29 Axial variation of the heat-transfer rate. Red = 8.3 x 10 4 ; a = 10°. (From Ref. 130.) 
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Fig. 31 Effect of angle of attack on 10° sharp cone heat transfer. (From Ref. 134; reprinted with 
permission of The American Institute of Aeronautics and Astronautics.) 
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Fig. 33 Comparison of heat-transfer coefficients. Moo = 7.4, a = 10°. (From Ref. 138.) 



